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ABSTRACT 

We report astrometric results of phase-referencing VLBI observations of 43 GHz SiO 
maser emission toward the red hypergiant VY Canis Majoris (VY CMa) using the Very 
Long Baseline Array (VLBA). We measured a trigonometric parallax of 0.83 it 0.08 mas, 
corresponding to a distance of 1.20lQ jg kpc. Compared to previous studies, the spatial 
distribution of SiO masers has changed dramatically, while its total extent remains 
similar. The internal motions of the maser spots are up to 1.4 mas yr~^, corresponding 
to 8 km s~^, and show a tendency for expansion. After modeling the expansion of maser 
spots, we derived an absolute proper motion for the central star of fix = —2.8 it 0.2 and 

= 2.6 lb 0.2 mas yr~^ eastward and northward, respectively. Based on the maser 
distribution from the VLBA observations, and the relative position between the radio 
photosphere and the SiO maser emission at 43 GHz from the complementary Very Large 
Array (VLA) observations, we estimate the absolute position of VY CMa at mean epoch 
2006.53 to be aj2ooo = 07'^22"^58?3259 ± 0?0007, 5j2ooo = -25°46'03^'063 ± O^'OIO. The 
position and proper motion of VY CMa from the VLBA observations differ significantly 
with values measured by the Hipparcos satellite. These discrepancies are most likely 
associated with inhomogeneities and dust scattering the optical light in the circumstellar 
envelope. The absolute proper motion measured with VLBA suggests that VY CMa 
may be drifting out of the giant molecular cloud to the east of it. 

Subject headings: astrometry — masers — parallaxes — proper motions — stars:individual 
(VY CMa) — stars:supergiants 
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INTRODUCTION 



VY Canis Majoris (VY CMa) is one of most massive and luminous red hypergiant stars in 
our Galaxy. At its "traditio nal" estimated dista nce of 1.5 kpc, based on the angular proxim- 
ity to the cluster NGC 2362 (jLada fc ReidI Il978l ). VY CMa wo uld be near the emp irical upper 
luminosity boundary in the Hertzsprung-Russel (H-R) diagram (jSchuster et al.ll2006l ). associated 
with high mass-loss and ejection phenomena. Because the red hypergiant phase represents a very 
short-lived evolutionary stage, with a timescale of only ~ 10^ years, the physical properties and 
evolutional state of this object have been investigated extensively from optical and infrared to radio 
wavelengths. 

VY CMa displays OH, H2O and SiO maser emission in its circumstellar envelope (CSE). Very 
Long Baseline Interferometry (VLBI) observations of SiO masers from oxygen-rich asymptotic gi- 
ant branch stars (e.g., Mira variables) have shown that the masers are situated only a few radii 
from the stellar surface, giving us insight into th e properties o f the inner CSE, which is important 
for the study of the evolution of late type stars (jBoboltal2005l ). Investigatior is of the distributions 



of the SiO masers show that the maser s typically form ring-like structures (IDiamond et al. 
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Creenhill et al. |[l995|;l Cotton et al.ll20o3), p ossibly with a rotational component (jBoboltz Sz Diamond! 
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VLBI observations of SiO maser emission from VY CMa in the J = 1^0, v = l and v = 2 
transitions showed them to be approximately spatially coincident, wit h extents of approximatel y 
80 mas on the sky and the strongest spots concentrated toward the east (iMiyoshi et al.lll994l . l2003l ) . 
VLBI maps of the f = l, 7 = 2—7-1 maser emission extend over about 1 00 x 80 mas in righ t 
ascension and declination, respectively, estimated at about 2—4 stellar radii (jShibata et al.ll2004l ). 



Knowledge of distance is very important to determine the physical properties of a star. Inferred 
stellar radii depend directly on distance and luminosities scale as the square of distance; moreover, 
distances indirect ly affect rnodeled stellar properties, such as effective temperature and surface 
gravity. Recently. IChoi et al.l (|2008al ) measured a trigonometric parallax for VY CMa of 0.88 it 0.08 



mas, corresponding to a distance of 1.14lQgg kpc, based on observations of H2O masers with the 
Japanese VLBI Exploration of Radio Astrometry (VERA) array. At this distance, VY CMa more 
comfortably falls below the theoretical maximum luminosity on the HR diagram. Our parallax 
measurements presented here are made with a different telescope, the Very Long Baseline Array 
(VLBA), observing a different maser species (SiO) and thus give an independent measurement of 
the distance and proper motion of VY CMa. 

While stellar positions for red giant stars have been measured with the Hipparcos satellite, 
these are generally far less accurate than for dwarf stars, since red giant stars are large, variable, 
often surrounded by copious dust and far away. Another way to loca te a red giant star is to observe 
circumstellar masers with radio interferometry (jBaudrv et al.lll984l ). Because masers usually sur- 
round the red giant star, the radio position of the central star can be estimated from the observed 
distribution and kinematics of the SiO masers, which can be linked to the positions of extragalac- 
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tic quasars by VLBI phase-referencing (jColomerill993l ). More directly, however, radio continuum 
emission from the stellar photosphere can be im aged by, for example, the Very La rge Array (VLA) 
using H2O or SiO masers as a phase-reference (|Reid &: Mentenlll99d . Il997l . 120071 ) . which provides 
a direct and accurate measurement of the maser distribution relative to the star. 

In this paper, we present the results from our multi-epoch radio interferometer observations 
of the SiO maser emission and radio continuum toward VY CMa. In ^ we describe the VLBA 
phase-referencing observations of SiO masers and the VLA observations of weak radio photospheric 
continuum and maser emission at 43 GHz made on the day after the second epoch of our VLBA 
observations. In ^ we compare the spatial distribution of SiO masers at different epochs and 
estimate their internal motions. In ^ we use the time variation of the positions of maser spots 
relative to a background source to determine a trigonometric parallax and absolute proper motion. 
In ^ we describe the procedures to determine the absolute position of the central star based 
on the quasi-ring-like structure of maser distributions from the VLBA observations and the VLA 
observations of radio photospheric emission using strong maser emission as phase-reference. In ^ 
after modeling the kinematics of masers, we derive an absolute proper motion of the star. In 
we compare the positions and proper motions with values measured by the Hipparcos satellite and 
VERA, and discuss the reasons for the difference and the origin of VY CMa. 



2. OBSERVATIONS 

2.1. VLBA Observations 

Our VLBI observations were conducted with the VLBA operated by the National Radio As- 
tronomy Observatory (NRAofl under program BR106. We observed the f = l, 7 = 1—7-0 SiO 
maser transition at a rest frequency of 43.12208 GHz toward VY CMa with 8-hour tracks on 2005 
October 20, 2006 April 16 and September 29, and 2007 April 14. This time sampling provides 
nearly maximum sensitivity for parallax detection and ensures that we can separate the secular 
proper motion from the sinusoidal parallax effect. We scheduled the observations so as to maxi- 
mize the right ascension (and not declination) parallax offsets for two reasons. First, for VY CMa 
the amplitude of the parallax signature in right ascension is greater than in declination. Second, 
since VY CMa is observed at low elevation angles, the uncertainty of declination meas urements 



would be expected to be considerably larger than for right ascension (jHonma et al.l 120081 ). 



We observed several extragalactic radio sources as potential background references for parallax 
solutions. The observing sequence was VY CMa, J0731-2341, VY CMa, J0720-2628, VY CMa, 
J0724-2515, VY CMa, J0725-2640 for the first epoch. From the first epoch observations, we 
found J0720— 2628 and J0724— 2515 too weak to be detected, and we dropped these sources from 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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the remaining observations. We switched between the maser target and background sources every 
40 s, typically achieving 30 s of on-source data. We used an SiO maser spot as the phase-reference 
source because it is considerable stronger than the background source and could be detected on 
individual baselines in the available on-source time. Table [1] lists the positions, intensities, source 
separations, LSR velocity of the reference maser feature and synthesized beam sizes. 



Table 1: Positions and Brightnesses 



Source 


R.A. (J2000) 


Dec. (J2000) 


Sp 


^sep 


P.A. 


Vlsr 


Beam 








(h m s) 


(° ' ") 


(Jy/beam) 


n 


n 


(km s^^) 


(mas 


mas 


°) 


VY CMa 


07 22 58.3283 


-25 46 03.075 


18 - 27 






33.3 


0.5 X 


0.2 (Q 


i -12 


J0725-2640... 


07 25 24.4130 


-26 40 32.680 


0.03 


1.1 


+ 34 




2.7 X 


0.6 C 


i -11 


J0731-2341... 


07 31 06.6680 


-23 41 47.869 


0.06 


2.8 


+136 




0.8 X 


0.3 C 


S +6 



Note. — The fourth and seventh columns give the peak brightnesses (Sp) and Vlsr of reference feature. The 
fifth and sixth columns give the separations (Osep) and position angles (P.A.) east of north between maser 
and background sources. The last column gives the FWHM size and P.A. of the Gaussian restoring beam. 
Calibrator J0725— 2640 is from the VLA pr ogram AR569, and it s position is corrected using information 
from th e VLBA Calibrator Survey (VCS) bv lXovalev et al.l |2007h . the information for J0731-2341 is from 
VCS bv lPetrov etaP (|2005h . 



We placed observations of two strong sources (J0730— 1141 and J0530+1131) near the begin- 
ning, middle, and end of the observations in order to monitor delay and electronic phase differences 
among the intermediate- frequency bands. The rapid-switching observations employed two adjacent 
bands of 8 MHz bandwidth and recorded both right and left circularly polarized signals. The two 
(dual-polarized) bands were centered at Local Standard of Rest velocities (Vlsr) of —33.6 and 22.0 
km for VY CMa. The SiO masers were contained in the second band. 

In order to do atmospher ic delay calibrati on, we placed "geodetic" blocks before and after our 
phase-reference observations (jReid et al.ll2009l ). These data were taken in left circular polarization 
with eight 8 MHz bands that spanned 480 MHz of bandwidth between 42.9 and 43.4 GHz; the bands 
were spaced in a "minimum redundancy configuration" to uniformly sample, as best as possible, 
all frequency differences. 

The data were correlated in two passes with the VLBA correlator in Socorro, NM. One pass 
generated 16 spectral channels for all the data and a second pass generated 256 spectral channels, 
but only for the single (dual-polarized) frequency band containing the maser signals, giving a 
velocity resolution of 0. 21 km s~^. The data calibration was performed with the NRAO AIPS 
package as described by iReid et al.l (|2009l ) . 
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2.2. VLA Observations 

Observation of tiie radio photosphere and maser emission toward VY CMa at 43 GHz were 
made with VLA in its largest (A) configuration under program AR595 on 2006 April 17. We used a 
dual intermediate frequency band setup with a narrow (6.25 MHz) band covering an LSR velocity 
range of —5 to 40 km s~^ for the f = 1,7=1—7-0 SiO maser line (rest frequency of 43122.08 MHz) 
and a broad (50 MHz) IF band centered ~ 50 MHz above the masers on a line-free portion of the 
spectrum. We also observed the SiO maser emission at high spectral resolution with several scans 
in spectral-line mode interspersed among the dual-band continuum observations. These covered 
the SiO line emission with a bandwidth of 6.25 MHz centered at 22 km s~^ using 128 spectral 
channels, which provided a channel spacing of 0.34 km s~^. 

A typical observing unit for our continuum observations consisted of a ~ 50 minute VY CMa 
scan, followed by a ~ 1.5 minute scan of, alternately, the quasar J0730-1141 and J0648-3044. 
For a spectral-line observation, we observed a calibrator J0730— 1141 for about 1.5 minutes and 
VY CMa for 6.5 minutes. Absolute flux density calibration was established by an observation of 
3C 286 (J1331-I-3030), as suming a flux density of 1.46 Jy at 43 GHz. Data calibration procedures 



are described in detail in lReid fc MentenI (jl997l . 120071 ). 



3. SPATIAL DISTRIBUTION AND INTERNAL MOTION OF MASER SPOTS 

Fig. [1] shows interferometer spectra of SiO maser emission observed with a long VLB A baseline 
at all four epochs. The SiO maser emission spans a I^SR range of about 6 to 42 km s~^. One can 
see that the flux densities of some maser features varied considerably over 1.5 yr. We selected a 
compact and relative stable strong maser spot at Vlsr of 33.3 km s~^ to serve as the phase-reference. 
The point-source response function (dirty beam) typically had a FWHM of 0.5 by 0.2 mas at a 
position angle of —12° east of north using uniform weighting. 
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Fig. 1. — Interferometer (scalar averaged cross-power amplitude over the full duration of the 
observation) spectra of the SiO masers toward VY CMa obtained with the long VLBA baseline 
from Pie Town (New Mexico) to Mauna Kea (Hawaii) at four epochs. The arrow points to the 
maser feature at Vlsr of 33.3 km s~^ which served as the phase-reference. 
(A color version of this figure is available in the online journal.) 
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Fig. 2. — Spatial distribution of SiO maser emission towards VY CMa observed at four epochs. 
Observation dates are indicated in the lower right corner of each panel. The Vlsr velocity of the 
maser spots is color-coded as indicated by the color bar on the right panel. 
(A color version of this figure is available in the online journal.) 
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Figure [2] shows the spatial distribution of SiO maser emission towards VY CMa relative to the 
reference maser spot at Vlsr= 33.3 km s~^. The total exte nt of the SiO maser spot distribution 



is about 80 mas, which is consistent with previous studies (jMivoshi et al.lll994l . |2003| ). while the 
spatial distribution changes in detail. We considered maser spots at different epochs as being from 
the same feature if their positions in the same spectral channel were coincident within 5 mas at 
different epochs, corresponding to a motion of less than 20 km s^^. Selecting only the brightest 
maser spot in each channel, we found nine features including at least two maser spots in adjacent 
channels that were detected at all four epochs (see Fig. [3]). 

The relative motions of features were calculated by performing a weighted least-squares fit for 
linear motion with respective to the reference feature. After removing the average relative motion of 
all maser features, the remaining (internal) motion vectors, indicated with arrows in Fig. [3l exhibit 
an expanding structure. These internal motions are also listed in Table [21 the maximal motion is 
measured for feature B. It is about 1.4 mas yr~^, corresponding to 8 km s~^ at a distance of 1.2 
kpc. 

Fig. S] shows the time variations of the positions of some maser spots relative to the reference 
spot. The masers in feature D exhibit approximate linear motions, while those in feature F show less 
regular motions. Indeed, in feature F, it is difficult to avoid possible misidentification of the same 
maser spots from one epoch to another, which introduces uncertainty in the kinematic parameters 
of those maser spots. Thus, for parallax measurements, only those with linear motions were used. 
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Fig. 3. — Averaged positions (circle and letter label) and internal motions (arrow) of the SiO maser 
features towards VY CMa. The average motion of aU features has been removed. The color bar 
denotes the VhSR range from 42 to 6 km s~^ of the maser features. The length and the direction 
of an arrow indicate the speed (given by the scale arrow in the lower left of the figure) and the 
direction of internal motion of a maser feature. The reference maser spot is in feature E, located 
at (0,0) mas, which is slightly different from the average position of maser spots in feature E. 
(A color version of this figure is available in the online journal.) 
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Table 2. Internal motions of maser features 



Region 


^LSR 












(km s~^) 


(mas yr 




(mas yr 




A 


15.0 


-0.187 ± 


0.029 


-0.511 ± 


0.003 


B 


18.5 


0.732 ± 


0.087 


1.359 ± 


0.007 


C 


23.7 


0.705 ± 


0.202 


0.621 ± 


0.055 


D 


8.5 


-0.664 ± 


0.007 


1.182 ± 


0.003 


E 


35.5 


0.240 ± 


0.013 


-0.377 ± 


0.011 


F 


22.0 


0.100 ± 


0.035 


-0.501 ± 


0.011 


G 


19.8 


0.418 ± 


0.007 


-0.455 ± 


0.012 


H 


39.8 


-0.735 ± 


0.070 


-0.573 ± 


0.057 


I 


27.6 


-0.610 ± 


0.114 


-0.745 ± 


0.143 



Note. — Motions fi^ ~ l^acosS and fiy ~ fis arc listed with the unweighted average motion removed. 
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Fig. 4. — Temporal changes of relative positions of the SiO maser spots in feature D and F. Circles, 
pentagons, squares and triangles denote the first to the last epoch, respectively. The error bars are 
10 times the formal position uncertainty. The Vlsr of each maser spot is indicated by the number 
close to the position at the first epoch. 

(A color version of this figure is available in the online journal.) 
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4. PARALLAX AND PROPER MOTION 

We fitted elliptical Gaussian brightness distributions to the images of strong maser spots and 
the extragalactic radio sources for all four epochs. The change in position of each maser spot 
relative to each background radio source was modeled by the parallax sinusoid in both coordi- 
nates (determined by a single parameter, the star's parallax) and a linear proper motion in each 
coordinate. 

Fig. [5] shows the maser reference channel images at all four epochs. One can see that the 
emission appears dominated by a single compact component, and there is no dramatic variation 
over the 1.5 year time span of our observations. The extragalactic source J0731— 2341 was relatively 
far (2.8°) from VY CMa, and it was only used to determine the absolute position of the maser 
reference spot. We only used J0725— 2640, separated by 1.1° from VY CMa, to determine the 
parallax. 

Fig. [6] provides images of background radio source J0725— 2640 using only the inner 5 VLBA 
antennas at all epochs. This source is dominated by a single component with peak brightnesses 
of 0.13, 0.04, 0.15 and 0.02 mJy beam~^ at our four epochs. The variation of the brightnesses is 
probably caused in part by significant flux loss owing to poor phase coherence in the two April 
observations, which were in evening hours when coherence is often poor. 
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Fig. 5. — Images of the reference SiO maser spot at Vlsr = 33.3 km s^^ in VY CMa. Observation 
dates are indicated in the upper left corner and the restoring beam (gray) is indicated in the lower 
left corner of each panel. Contour levels are spaced linearly at 3.0 Jy beam^^. 
(A color version of this figure is available in the online journal.) 
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Fig. 6. — Images of extragalactic radio source J0725— 2640 used for the parallax measurements of 
VY CMa. Observation dates are in the upper left corner and the restoring beam (gray) is in the 
lower left corner of each panel. Contour levels are spaced linearly at 25, 25, 10 and 5 mJy beam~^ 
for the four epochs in chronological order. 

(A color version of this figure is available in the online journal.) 
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As mentioned in ^ the apparent motions of the maser spots can be comphcated by a combi- 
nation of spectral blending and changes in intensity. Thus, for parallax fitting, one needs to find 
stable, unblended spots and/or use many maser spots. Thus, we first fitted a parallax and proper 
motion to the position offsets for each maser spot separately (with respect to the background ra- 
dio source J0725— 2640). In Fig. \7\ we plot the position of the reference maser spot relative to 
the background radio source as an example, with superposed curves representing the model maser 
tracks across the sky. This is one of the most compact and unblended spots that was detectable at 
all four epochs. 




3 2 1 -1 -2 -3 2006 2007 2006 2 

East OHset (mas) Epoch (year) Epoch (year) 



Fig. 7. — Parallax and proper motion data (circles) and best-fitting models (lines) for the reference 
maser spot at Vlsr of 33.3 km s~^. Plotted are positions of the maser spot relative to the extra- 
galactic radio source J0725— 2640. Left Panel: Positions on the sky for the first and last epochs 
are labeled. The expected positions from the parallax and proper motion fit are indicated with 
(crosses). Middle Panel: Eastward (solid lines) and northward (dashed lines) offsets and best-fitting 
models versus time. Data for the eastward and northward positions are offset vertically for clarity. 
Right Panel: Same as the middle panel, except the best-fitting proper motion has been removed, 
displaying only the parallax signature. 

(A color version of this figure is available in the online journal.) 
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Fig. [8] shows the estimated parallaxes and proper motions for all measured maser spots. While 
the parallaxes should be identical (within measurement uncertainties), the proper motions are 
expected to vary among the spots owing to internal motions of ~ 10 km (1.8 mas yr~^ at 
a distance of 1.2 kpc). While most of the maser spot parallaxes show good internal consistency, 
the dispersion of parallax estimates over the entire ensemble is considerably larger than the formal 
errors would suggest. This is caused by residual systematic errors affecting the fits, which originate 
in the complexity and evolution of blended spectral and spatial structure for some of the masers. 

We discarded all parallax solutions that had formal uncertainties large than 150 /xas, corre- 
sponding to relative uncertainties of 18% (for a parallax of 0.83 mas). Some parallaxes from region 
G appear to be outliers, and we discarded the results from this region. The remaining fits generally 
yield internally consistent parallaxes. Most of the parallax results are distributed in a range defined 
by n lb 2(Tn, where 11 and an are the average and standard deviation of parallaxes, respectively. 
Table [3] lists the results of the parallax and proper motion fits for all remaining maser spots. 

Since one expects the same parallax for all maser spots, we did a combined solution (fitting 
with a single parallax parameter for all maser spots, but allowing for different proper motions for 
each maser spot) using all the remaining maser spots (indicated with squares in the top panel 
of Fig. [8j). The combined parallax estimate is 0.830 it 0.079 mas, corresponding to a distance of 



1.20"'"q ^3 kpc, which is consistent with the VERA result for the H2O masers in VY CMa (iChoi et al 



2008al ). The quoted uncertainty is the formal error multiplied by ^/n (where n is the number 
of maser spots used in the final parallax fit) to allow for the possibility of correlated position 
variations for all the maser spots. This could result from small variations in the background 
source or fr om unmodeled a tmospheric delays, both of which would affect the maser spots nearly 



identically (jReid et al.ll2009l ). The average proper motion of all selected maser spots is —2.21 ±0.06 
mas yr~^ eastward and 2.29 ± 0.30 mas yr~^ northward, which is reasonably compatible with the 
VERA results (—2.21 it 0.16 mas yr~^ eastward and 1.02 ± 0.61 mas yr~^ northward). 
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Fig. 8. — Individual solutions for parallax and proper motion of all measured maser spots. Regions 
of maser spots are indicated with letters near the bottom of the top panel. The dash-dotted lines 
indicate the means and the dotted lines indicate the mean ±2 standard deviations for the parallaxes. 
(A color version of this figure is available in the online journal.) 
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Table 3. Parallax and proper motion fits 



Region 


Ch. 


^LSR 




Parallax 






fJ.y 






(km s"-*-) 




(mas) 


(mas 






(mas 


y ) 


A 
A 


iDl 


1 c; n 
10. U 


U 




-1- n ci'\R 


-2.625 


± 


U.UDo 


2.196 


± 0.365 




lOZ 




1 


.007 


i 0.044 


-2.323 


it 


U.U t o 


2.252 


it 0.457 


L) 




o.O 


U 


■ O-L 1 




-3.081 


it 


U.zzo 


3.972 


it 0.466 




1 no 


Q Q 

o.O 


U 




-1- 1 


-3.088 


it 


U.zDU 


3.982 


it 0.454 


III 


oo 


oo.O 


n 
U 


.770 


i 0.129 


-2.437 


it 


U.Zol 


2.164 


it 0.660 




OO 


QQ 1 
OO. 1 


U 




-1- n 1 "^7 


-2.544 


it 


U.z40 


2.385 


it 0.455 




Oo 


o / .y 


U 


.763 


it 0.113 


-2.492 


it 


V.ZKJZ 


2.346 


it 0.506 




/ 


O / .D 


U 


.768 


lb 0.087 


-2.460 


it 


U.lOO 


2.325 


it 0.553 




Oo 


o / .4 


U 


777 
• III 


-t- n HQ'? 


-2.475 


it 


U.IDO 


2.310 


it 0.538 




Oy 


O / .Z 


n 
U 


.793 


it 0.143 


-2.431 


it 


U.ZO / 


2.296 


it 0.517 




DU 


o / .U 


U 


. ouu 


-t- n 1 '^7 
m u. lo I 


-2.417 


it 


U.Z40 


2.296 


it 0.526 






Q^; 1 
oo. 1 


U 


.779 


i 0.082 


-2.406 


it 


U.140 


2.260 


it 0.453 




DO 


oo.y 


U 


.824 


it 0.053 


-2.334 


it 


U.uy4 


2.242 


it 0.459 




DO 


oO. ( 


U 


R97 


m u.u / u 


-2.173 


it 


U.lOD 


2.217 


it 0.470 




D / 


oO.O 


U 


.812 


it 0.080 


-2.107 


it 


U.140 


2.201 


it 0.456 




Do 


oO.o 


U 


.770 


it 0.079 


-2.142 


it 


U. 141 


2.218 


it 0.463 




Dy 


oO.U 


U 


.786 


it 0.042 


-2.100 


it 


U.U ID 


2.225 


it 0.448 




f U 


o4.o 


U 


.829 


it 0.036 


-2.001 


it 


U.UD4 


2.222 


it 0.503 




f 1 


o4.D 


U 


.844 


it 0.012 


-2.025 


it 


U.Uzz 


2.216 


it 0.576 




72 


34.4 


0, 


.840 


it 0.078 


-2.051 


it 


0.139 


2.262 


it 0.610 




73 


34.2 


0, 


.802 


it 0.070 


-2.080 


it 


0.124 


2.321 


it 0.585 




74 


33.9 





.783 


it 0.057 


-2.071 


it 


0.101 


2.350 


it 0.540 




75 


33.7 


0, 


.798 


it 0.076 


-2.059 


it 


0.135 


2.366 


it 0.515 




76 


33.5 


0, 


.855 


it 0.050 


-2.094 


it 


0.089 


2.386 


it 0.488 




77 


33.3 


0, 


.857 


it 0.028 


-2.123 


it 


0.050 


2.411 


it 0.475 




78 


33.1 


0, 


.804 


it 0.045 


-2.240 


it 


0.080 


2.442 


it 0.452 




79 


32.9 


0, 


.745 


it 0.074 


-2.280 


it 


0.133 


2.462 


it 0.388 




80 


32.6 


0, 


.759 


it 0.113 


-2.175 


it 


0.203 


2.481 


it 0.373 




81 


32.4 


0, 


.888 


it 0.078 


-1.792 


it 


0.140 


2.523 


it 0.409 




82 


32.2 


0, 


.983 


it 0.038 


-1.453 


it 


0.068 


2.557 


it 0.436 




83 


32.0 


0, 


.928 


it 0.034 


-1.485 


it 


0.061 


2.529 


it 0.369 


F 


129 


22.0 


0, 


.834 


it 0.057 


-2.322 


it 


0.101 


2.236 


it 0.623 




130 


21.8 


0, 


.903 


it 0.027 


-2.252 


it 


0.048 


2.322 


it 0.530 




131 


21.6 


1, 


.021 


it 0.036 


-2.102 


it 


0.064 


2.261 


it 0.568 




132 


21.3 


1, 


.073 


it 0.009 


-1.993 


it 


0.016 


2.237 


it 0.598 


I 


103 


27.6 


0, 


.757 


it 0.125 


-2.979 


it 


0.224 


2.208 


it 0.649 




104 


27.4 


0, 


.784 


it 0.079 


-2.997 


it 


0.142 


2.240 


it 0.794 


Combined 






0, 


.830 


it 0.079 


-2.210 


it 


0.060 


2.290 


it 0.300 



Note. — Absolute proper motions arc defined as fi^ — fJ-acosS and fiy — fJ-s- 
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5. ABSOLUTE POSITION OF CENTRAL STAR 
5.1. Absolute position of VY CMa derived from VLB A observation 



Based on the circumstellar maser distribution, it is possible to locate the central star with 
reasonable accuracy. For example, for a ring-like distribution of maser emission, the star should be 
located at the center of the ring. Th at this is indeed the c ase has been obs e rvatio nally verified for 
AGE stars by iReid k MentenI (jl99d ) for H2O masers and iReid fc MentenI (120071 ) for SiO masers. 
Previous VLBI observatio ns of SiO masers towards VY CMa have shown that the masers appear 
clumpy and asymmetrical (jMiyoshi et al.lll994l . 120031 ). Our observations (see Fig. [2]) show that the 
SiO masers toward VY CMa exhibit only a partial ring-like structure. 

There are "spoke-like" (linearly distributed) maser features at all four epochs; one of the 
best examples is found in the group of features which include the reference spot (see Fig. [9]). These 
maser features, appear to be composed of "spokes" of gas flowing outward from the c entral star and 



mostly have decreasing radial velocities with increasing the distance from the star (lYi et al 



2005 



Matsumoto et al.ll2008l ). There are several spoke-like maser features in VY CMa and they probably 
"point back" to the star. From these spoke-like features and the ring-like maser distribution, we 
estimate the position of the central star relative to the reference spot to be (Ax = —28, Ay = +10 
mas) in east and north directions with an uncertainty of about 10 mas. Based on the combined 
parallax fit in ^ we determine the position offset (Ax = —4.8 it 0.1, Ay = 1.6 it 0.2 mas) of the 
reference spot at mean epoch 2006.53 of our observation. Combining these offsets, we estimate the 
absolute position (aj2ooo = 07'^22™58?3259, 5j2ooo = -25°46'03^'063) of the central star with an 
uncertainty of about 10 mas. 
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Fig. 9. — Positions of SiO masers for 2006 April 16 and inferred position of t he central star (crosg ) 
with a radius (circle) of 10 mas estimated from mid-infrared observations (iMonnier et alJl2000l ). 
Each maser spot is represented by a filled circle whose area is proportional to the logarithm of the 
flux density. The dotted lines labeled with a number from "spoke-like" maser features point back 
to the central star. The Vi,sR velocity of the maser spots is color-coded as indicated by the color 
bar on the right side. 

(A color version of this figure is available in the online journal.) 



5.2. SiO masers and radio photosphere from VLA observation at 43 GHz 



From mul ti-wavelength VLA observat ions of a sample of long period Mira and semi-regular 
variable stars, iReid &: MentenI (Il997l . 120071 ) detected emission from their radio photospheres. The 
approach of using a strong narrowband signal (maser emission) in one observing band to calibrate 
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the phase (and amphtude if necessary) of a weak br oadband signal (radio c ontin uum emission from 
the radio photosphere) was described in detail by iReid &: MentenI (|l997l . 120071 ) . In the fohowing 
we refer to the two datasets as narrow band (NBD) and broad band (BBD) data. 

In order to register the SiO masers relative to the radio photosphere, we first measured the 
position of the radio photosphere from the BBD relative to the maser emission in the (pseudo- 
continuum) NBD. Next, we aligned the maser emission from the spectral- line data by producing 
a map of "simulated" NBD (hereafter SNBD) from the line data. Because the NBD and SNBD 
cover the same velocity range, we can align the maps from the NBD and SNBD by comparison of 
positions of peaks in each map. This allows the positions of the emission in individual channels to 
be registered to the NBD, and then to the BBD (radio photosphere) emission. 



Following the procedures described in detail in iReid &: MentenI (119971 ) . once the data were 
"cross-self-calibrated", both NBD and BBD were imaged with the AIPS task IMAGR. We per- 
formed Gaussian position fits to the peak in the BBD and NBD, using AIPS task JMFIT, and 
obtained the relative position (Ax = —16 it 2, Ay = —6 it 2 mas). 

The line data were self-calibrated by choosing a channel with strong emission as a reference, 
and the resulting phase and amplitude corrections were applied to the other channels. From the line 
data we produced a spectral-line data cube and then we produced a NBD-like map by averaging 
all the channels in the cube. We fitted the two strongest peaks in the NBD and SNBD maps and 
calculated a relative position of (Ax = —8, Ay = +24 mas) with an uncertainty less than 1 mas. 
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Fig. 10. — SiO maser distribution around VY CMa. All offsets are relative to the best fit position 
of the peak in the map of the strongest maser emission at Vlsr of 22.1 km s~^. All symbols 
are overlaid on contours of the 20 mas resolution image of the integrated narrowband emission 
(SNBD). Contour levels are spaced linearly at 5 times the rms noise level of 0.69 Jy beam~^. All 
the maser spots are marked with circles^ whose sizes are proportional to the logarithm of their flux 
densities. The best fit position of the radio photosphere (BBD) is marked as a cross. The dotte d 
circle presents the boundary of the central star with a radius of 10 mas from iMonnier et al.l (j200d ) . 
The Vi^sR of the maser spots is color-coded as indicated by the color bar to the right. 
(A color version of this figure is available in the online journal.) 



Finally, we estimated a position offset (Ax = —24, Ay = +18 mas) of the weak continuum 
source (radio photosphere) relative to the maser spots (reference spot) with a uncertainty of 3 
mas. Thus, all the maser spots were registered to the radio photosphere at a similar accuracy (see 
Fig. do]). 

However, to estimate absolute position of the continuum source, we have to align the maser 
spots with those from the VLBA observations. Due to the significantly different angular resolution 
and sensitivity of the VLA and VLBA, there is no guarantee that the two images should be identical. 
Note that the VLA has a 50 mas synthesized beam, which for any channel essentially gives a centroid 
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position. If a single spectral channel contains emission from across the source, the VLA position 
would be falsely interpreted as near the center of the distribution. 
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Fig. 11. — Comparison of SiO maser distributions from the VLBA (circle) and VLA (square), 
after cross registration (see text). The observation epochs are also labeled. The position of the 
star determined from the VLA radio photosphere observations is indicated by the "+" and from 
the VLBA "spoke- like" maser spot distribution by the "x" symbols. 
(A color version of this figure is available in the online journal.) 



Nevertheless, the SiO maser distributions from VLA and VLBA data bear significant resem- 
blance. In order to align these two maps, we assumed that the centroid positions of the strongest 
maser spots are consistent. For the VLA map, the strongest maser spot was chosen as the phase 
reference and, hence, is at the map origin (see Fig. [TOj) . For the VLBA map, the phase reference 
was not the strongest spot, and the reference spot is at (-|-3, —10) mas (see Fig. [9]). After shifting 
all positions from the VLA map by this offset, we find reasonable agreement with the VLBA and 
VLA maser spot maps. Since we obtained an offset of (—7, +2) mas for the central star (imaged 
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with the VLA) relative to the VLBA spots maps, we now can locate the star on the VLBA spot 
maps to about 10 mas uncertainty. This is shown in Fig. [TTJ 



6. THE KINEMATIC MODEL 



6.1. Expansion 

Combining the proper motion (/i^., fXy) and parallax with the of each maser spot, we can 



derive its 3-dimensional velocity vector Vi{vix,Viy,Viz) and its variance 



4. 



turb 



(1) 

(2) 
(3) 



where k is used to convert proper motion to velocity adopting the distance of 1.2 kpc from ^ and 
vturb is a "turbulent" velocity component, as s umed to be ~ 3 km s~^. The Vlsk of the star (Vlsr*) 
is adopted as 22 km from iMenten et al.l (j2006l ). 



The model velocity vector for an expanding flow is 



exp I 



(4) 



where, v^xp is an expansion speed modeled as 



^exp 



(5) 



Vo{vox,voy,voz) is the systemic velocity (presumably of the exciting star) and ri{Ri — Rq) is the 
position vector of a maser spot Ri{xi,yi, Zi) relative to the expansion center RQ{xo,yo, zq). The 
sky positions {xi,yi) of the maser spots are relative to the reference maser spot which is located at 
(x = 0,y = 0), and the radial position z is relative to the expansion center assuming z = 0. The 
sky positions {x,y) are converted from mas to AU adopting a distance of 1.2 kpc. 

The three components of velocity for all maser features in ^are taken as observables, and the 
weight of each observable is assig ned the reciproca l of its variance. The fitting procedur e is similar 



to that used in previous studies (jReid et al. 



1988 



Gwinn et al 



1992 



Imai et al.l 120001 ): however 



unlike most other expanding maser models, we know the distance and vqz and do not have to 
solve for them. Apart from the global parameters listed in Table [H the Zi of each maser spot 
is also solved for as a free parameter. In order to restrict the parameter estimates to physically 
reasonable values as well as to incorporate our knowledge of parameters xq and yo, we added a priori 
information for these two parameters by adding extra data-like equations, one for each parameter, 
to the observables, and weighting the "pseudo-observables" with their estimated uncertainties. This 
allowed us to incorporate the information that {xq = —33 it 12, yo = 12 it 12 AU) from ^ 
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The global parameters and the Zi of each maser spot were adjusted to minimize the sum of 
the squares of the weighted residuals using the Levenberg-Marquardt method. We first performed 
the model fitting using ah the maser spots, and then discarded those maser spots (6 of 60) with 
large residuals (> Scr), and estimated the parameters again; the best fitted global parameters are 
listed in Table HI As Fig. [12] shows, the expansion velocities are much smaller than the systemic 
velocity, indicating that the absolute proper motions of the maser spots are dominated by the 
systemic velocity of —16 it 1 km s~^ eastward and 15 it 1 km s~^ northward (corresponding to 
= —2.8 lb 0.2 mas yr~^ and /Xy = 2.6 it 0.2 mas yr~^), which is consistent with the proper motion 
estimated in 21 within 3cr. 



Table 4: Kinematic Model of VY CMa SiO masers 



Parameter 


Value 


Uncertainty 


Units 




Comments 


Xq 


-35 


± 5 


(AU) 




Center of expansion 


yo 


13 


± 5 


(AU) 






Vox 


-16 


± 1 


(km s~ 




Systemic velocity 


VOy 


15 


± 1 


(km s~ 






VexpO 


2 


± 1 


(km s~ 




Expansion speed at the origin 


a 


0.09 


± 0.01 


(km s" 


-1 AU-i) 


Expansion speed gradient 
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Fig. 12. — Positions and velocities for maser spots. The maser spots are located at the apexes of 
the red cones. The red cone lengths and orientations show the speed and direction of their velocities 
after removal of the systemic velocity Vq, which is indicated with blue cone. The expansion center 
is located at the apex of the blue cone with the indicated error bar in X and Y directions. 
(A color version of this figure is available in the online journal.) 



6.2. Radial Spoke-like Maser Features 



> 1 




-10 -5 5 10 15 20 

Fig. 13. — Angular distances versus radial velocities (Vlsr~" ^lsr*) of spoke-like maser features. 
Different markers denotes different maser features from Fig. O Feature 4 is divided into three parts 
(a, b and c). 

(A color version of this figure is available in the online journal.) 

As mentioned in ^5.11 there are several spoke-like maser features in VY CMa. Fig.ll3lshows the 
distance-velocity structure of six of these features, labeled as in Fig. [9l These spoke-like features 
have clear velocity gradients, and most of them have radial velocities (Vlsr — ^SR*) > and 
decreasing Vlsr with increasing angular distance from the center. This could happen in two ways: 
(1) the spoke-like features have motions outward from the star (i.e., outflow) on the far side and 
are decelerating; (2) the spoke-like features have motions inward toward the star (i.e., infall) on 
the near side and are accelerating. Taken into consideration the motions shown in Fig. 3, we believe 
most of the spoke-like features are on the far side and are decelerating. 



In order to explore the natu re of the radial ve l ocities of the observed spoke-like features, we 
adopt the ballistic-orbit model of iMatsumoto et al.l (|2008l ) to fit the distance- velocity structure of 
the spoke-like feature. This model assumes that the maser spots are moving in a radial direction 
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with acceleration or deceleration due to the stellar gravity. 



The relation between the projected distance from the center of the star, r, and the line-of-sight 
velocity, v, is written as 

ro cos i 



1 + 



(6) 



2GM \smi. 

where v = Vlsr — Vlsr*, G is the gravitational constant, M the stellar mass, i (7^ 0) the incli- 
nation angle of the maser feature projected onto the sky plane, and ro is the distance of apoc- 
en ter (where f = 0) from the center of the star. However, unlike the ballistic-orbit mo del fitting 



m 



Matsumoto et al.l ( 20081 ) . we adopt a stellar mass of 25 Mq from lMonnier et al.l ( 200d ) . while we 



think ro is unknown and may be different for each spoke-like feature. Thus, we solve for tq, which 
is converted to angular units adopting a distance of 1.2 kpc, and i for each spoke-like feature. 



> 


rO= 


70.7 +/- 


1.4 (mas) i= 


54.7 +/- 


0.8 (deg) 


□ 


rO= 


35.5 +/- 


0.2 (mas) i= 


12.0 +/- 


0.8 (deg) 





rO= 


65.4 +/- 


1.3 (mas) i= 


45.6 +/- 


1.8 (deg) 





rO= 


46.7 +/- 


1.5 (mas) i= 


54.7 +/- 


1.3 (deg) 


V 


rO= 


35.7 +/- 


1.1 (mas) i= 


54.7 +/- 


1.3 (deg) 




rO= 


153.7+/- 


- 6.7 (mas) i= 


= 83.0 +/- 


- 0.3 (deg) 




LSR 



LSR 



10 

(km/s) 



Fig. 14. — Fitted ballistic orbit (dashed lines) for the spoke-like feature (markers) in the distance- 
velocity diagram. Different markers denote different maser features labeled in Fig. [131 The fitted 
parameters and their uncertainties are labeled in the top right corner. 
(A color version of this figure is available in the online journal.) 
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Fig. [T3] shows the observed and modeled distance- velocity structures of all the spoke-like fea- 
tures (except 6 and 4b which do not follow a ballistic model) using Eq. [6l The distance-velocity 
structure should be symmetric about v = 0, which means blue shifted {v < 0) maser feature should 
increase their li,sR with increasing distance. Since most of the spoke-like features can be reasonably 
well modeled with ballistic orbits, this suggests that most of the spoke- like features "point back" 
to the star. 
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DISCUSSION 



Images of some AGB stars such as IK Tauri demonstrate that SiO masers may form circular or 
ehiptical rings of emission. Also the line-of-sight velocity structure of some SiO masers can have an 
apparent axis of symmetry consiste nt with the elo ngation axis of the maser distribution, suggestive 
of rotation of the SiO maser shell (lBoboltal2005l . and references therein). However, for VY CMa 
such symmetry is not evident in the spatial or vel ocity distribution of th e mas er s. Large scale 



optica l and near- inf r ared ir nages of VY CMa (e.g. iKastner &: Weintraubl (jl998l ): iMonnier et al 



(jl999l ): ISmith et all (120011 ): Humphreys et al.l (120071 )) display a possible northeast-southwest axis 
of symmetry, which might suggest a possible bipolar structure. However, these observations are 
from dust emission much further from the star than the circumstellar SiO masers. After subtracting 
the average internal motion, the SiO masers show a tendency for expansion from th e central star, 
which is similar to the 22 GHz H2O masers measured by VERA (jChoi et al.ll2008bl ). Observation 
of SiO {v = 0, J = 1 —> 0) maser emission towards VY CMa with VLA shows a bipolar outflow 
along the line of sight (IShinnaga et al.l 12004 ). consistent with the 3-dimensional kinematics of the 
H2O masers from lChoi et al.l (|2008bl ). However, there is no strong evidence for bipolar outflow but 
a slow quasi-spherical outflow displayed in Fig. [12] from our VLBA observations. 

The 3-dimensional expansion model implies that the radial velocities of the maser features 
will increase with increasing distance from the star. If the maser features are radially aligned as 
the ballistic orbit model assumes, then the expansion model seems counter to the ballistic model. 
However, the ballistic-orbit model provides only one of the possible explanations for the spoke- 
like features. As mentioned in ^6.2|, there are still spoke-like features that cannot be modeled 
with a ballistic orbit. A more realistic model might involve a hydrostatic inner envelope and /o r 
acceleration driven by pulsation or giant convective cells first introduced by ISchwarzschildl (| 19751 ). 
For the expansion model, a very small value of the expansion speed gradient (see Table |4]) is rather 
weak evidence for acceleration in outward motion. Furthermore, even though we know the distance 
of the star, based on the combined parallax of the maser spots, the distances of the maser features 
from the star are unclear (i.e., coordinate component along the line of sight, z) and can only be 
solved for as parameters, which are dependent on the model. Nevertheless, the uncertain z's of 
the maser features for the expansion model will not significantly affect the derivation of the proper 
motion of the central star, which is one of the most important results in this paper. 

Our accurate parallax based on SiO masers measured with VLBA agrees with that determined 
from H2O masers with VERA (see Table [5]). The small difference in proper motions of the center of 
expansions for the two maser species is not strongly statistically significant and, even if real, might 
be caused by small errors in determining the center of expansion for the different maser species 
which are located in different places in the circumstellar envelope. 



Lada fc Reidl (j 19781 ) suggested that VY CMa and the open cluster NGC 2362 are at nearly the 



same distance, which was estimated to be 1.5 kpc from main-sequence fitting of NGC 2362. Their 
argument rested on the fact that VY CMa was located at the apex of an arc of bright emission that 
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is visible in red-light images and had a I^sr similar to that of the molecular cloud complex and 
the stars of the cluster NGC 2362. Lada & Reid also pointed out that there were other possible 
ionizing sources for the arc, including r CMa (HD 57061); UW CMa (HD 57060, 29 CMa); and 
NN CMa (HD 58011). These stars were estimated to be 1.1 kpc distant, based on photometric 
measurements, and significantly closer than NGC 2362. Since our trigonometric parallax is closer 
to 1.1 than 1.5 kpc, this suggests that the main-sequence fitting at the time was less accurate than 
the photometric distances. 

The main-sequence fitting distances are nearly 40 years old and one might worry that bet- 
ter calibrations would resol ve the discrepancy. Howeve r, more rec ent estimates of th e distance 
to NGC 2362 are 1.6 2 kpc dMermilliod Maederl llQsd'l. 1 .42 kpc (|Brown et aLlllOsd l. 1.48 kpc 



(jMoitinho et ahlboOll ) and 1.36 kpc ( Mavne fc Navlorl I2OO8I ). These are all reasonably consistent 



with the 1.5 kpc value previously used and still different from the parallax values determined in 
this paper and by VERA. For an individual star, a photometric distance is highly deperi d ent o n 
the spectral type and extin ction model. Compared to 1.1 kpc quoted by iLada &: ReidI (jl978l ). 
Kaltcheva fc HilditchI (j2000l ) estimate a significant different photometric distance of 2.86 kpc for 
NN CMa and 1.51 kpc for UW CMa. These are significantly larger than 1.2 kpc distance to VY 
CMa reported in this paper. Provided that VY CMa is indeed associated with NGC 2362 and NN 
CMa or UV CMa, this suggest that systematic uncertainties for the distance estimated using either 
main-sequence fitting or the photometric method is at least ~ 20%. 

The absolute position of the central star from our VLBA and VLA measurement at 43 GHz 
is consistent with that from VERA. However, the position and proper motion from Hipparcos are 
significantly different with the values determined from radio measurements (see Table [5|) . Using 
the position and proper motion from radio measurement, we predicted the position of VY CMa 
at the same epoch as the Hipparcos catalogue, then compared the position derived from optical 
and radio observations, and found a large discrepancy of tens of ma s. VY CMa is embedded in 



an asymmetric du s t reflection nebula, and recent images from optical (IKastner &: Weintraub 



1998 



Smith et al.l I2OOII : [Humphreys et al.l 120071 ) and near-IR observations (iMonnier et al.l 1 19991 ) show 
that the dust shell is time variable. Since most of the optical radiation coming from the star is 
reprocessed by the surrounding dust, a likely explanation for the difference between Hipparcos and 
the radio results is that the observed optical light is mostly scattered by circumstellar dust and a 
reliable stellar proper motion determination at optical wavelengths is practically impossible. 

Another possibility for a difference between the VLBA and Hipparcos proper motions would 
be that an inhomogeneity on VY CMa's surface contributes to the Hipparcos proper motion. For 
example, the early M-type supergiant a Orionis' ultraviolet appearance, as resolved w ith the Hubble 
space telescope, shows a conspicuous, compact "hot spot" (ICilliland fc: Dupreelll996l ). with a > 200 
K higher temperature than th e rest of the surface. One of their interpretations for this spot is in 
the context of a conjecture by ISchwarzschildl (119751 ) that a red giant star's surface only may show 
a few hundred convective cells, with large temperature diffe rences (of up to 1 00 K) between the 
hot and cool portions of a convective element. In addition, ISoker &: ClavtonI (|l999l ) investigated 
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the structure of cool magnetic spots in the photospheres of evolved stars and found the spots wil l 
cause the AGB star to appear asymmetrical. Adopting a diameter of 20 mas ([Monnier et al.1120001 ) 
for VY CMa, could the apparent proper motion from Hipparcos be the result of a hot spot rotating 
with the star's surface? Unfortunately, the rotation velocities of red supergiant are expected to 
be of order 1 km s~^, much lower the measured proper motion discrepancy. Could the random 
appearance and disappearance of different independent spots mimic the Hipparcos proper motion? 
Given that Hipparcos' proper motion determinations are typically based on 30 measurements within 
the mission lifetime of 3 yr, this is conceivable. 

Fig. [15] shows the geometry of VY CMa's optical environment overlaid with radio emission con- 
tours measured with the Atacama Pathfinder Experiment (APEX) telescope bv lSchuller &: Menten 
(|201ll ). The V-shaped surrounding material highlighted as Ha emission, which is either external to 
VY CMa (i.e., interstellar material) or comes from VY CMa (stellar mass loss), suggests relative 
motion of gas and the star. Taking our VLBA absolute proper motion, one explanation could be 
that VY CMa is drifting out of the molecular cloud to the east of it. 

Assuming that VY CMa was formed in the region where current signs of star formation are 
found (i.e., to the south-east of the molecular cloud as indicated by the color contours in Fig. [T5|) 
and that the absolute proper motion of the molecular cloud is small, then the travel time of VY 
CMa to its current location is about 0.5 Myr. This might be reasonable since VY C Ma is a very 
massiv e, and hence short-lived, star. However, adopt ing a luminosity of 3 x 10^ Lq from lchoi et al.l 
( 2008a ) and an effective temperature of 3650 K from lMassey et aL ( 2006 ). the location of VY CMa 
on an HR diagram is consistent with an age of of 8.2 Myr, based on an evolutionary track of a 
25 Mf7) star of solar me tallicity with an initial rotational velocity of 300 km s~^ at its equator 



(jMevnet fc Maederll2003l ). This age discrepancy between the motion time and age of VY CMa may 
imply that the relative proper motion between the molecular cloud and VY CMa is much smaller 
than its absolute proper motion or that the stellar evolution model is not appropriate for VY CMa. 



Table 5: Absolute position and proper motion of VY CMa from different telescopes 



Telescope 


Epoch 




R.A. (J2000) 




Dec. (J2000) 


fix 
















(h m 








(° ' ") 




(mas 






(mas 






Hipparcos" 


1991.25 


07 


22 


58.3251 


± 


0.0001 


-25 


46 03.180 ± 


0.003 


+8.86 


± 


1.34 


+0.75 


± 


3.25 


Hipparcos'' 


1991.25 


07 


22 


58.3251 


± 


0.0002 


-25 


46 03.176 ± 


0.003 


+5.72 


± 


2.01 


-6.75 


± 


4.47 


VERA'= 


2006.82 


07 


22 


58.3264 


± 


0.000? 


-25 


46 03.066 ± 


0.00? 


-2.09 


± 


0.16 


+ 1.02 


± 


0.61 


VLBA'' 


2006.53 


07 


22 


58.3259 


± 


0.0007 


-25 


46 03.070 ± 


0.010 


-2.21 


± 


0.06 


+2.29 


± 


0.30 



Perrvman et all (119971) . 



.van LeeuwenI (120071) . 

Choi et al.l (|2008allbh . where the position uncertainty is unknown. 



This paper. 




Fig. 15. — Ha emission {gray scale) superposed on the radio continuum emission at 870 /im with 
the Large Apex BOlometer CAmera (LABOCA) (green contours), and the CO(3-2) emission hne 
is integrated over 15-21 km s~^ {red contours) and over 21-27 km {blue contours). VY CMa is 
located at the bright rim of the V-shaped structure of the cloud complex and also a point source 
in the continuum emission contours. The arrow originating at VY CMa represents a motion over 
0.5 Myr adopting the absolute proper motion derived in this paper. The motion uncertainty is 
indicated as an error ellipse with a confidential level of 95% centered on the head of the arrow; 
the dashed error ellipse denotes the probable original position of VY CMa 0.5 Myr ago. Three 
candidates fo r the sources of ionization in this region: Tau CMa, UW CMa and NN CMa from 



Lada fc ReidI (j 19781 ) are labeled. The apparent size of ~ 8' of the nearby 0-star cluster NGC 2362 



centered on Tau CMa is indicated as a circle. 

(A color version of this figure is available in the online journal. 
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8. CONCLUSIONS 

We have measured the trigonometric parallax and proper motion of VY CMa from VLBA 
observations of the variations of relative positions between 43 GHz SiO masers and the background 
source J0725— 2640. The parallax of 0.83 it 0.08 mas provides a distance of 1.20jjn kpc, whic h 



is consistent with that measured with 22 GHz H2O maser using the VERA by I Choi et al.l (l2008al ). 
There can be little doubt that VY CMa is nearer than the value of 1.5 kpc generally used in the 
literature. 

Using the VLA detection of the radio photosphere, using the SiO masers as a phase reference, 
we determined the position of the central star relative to the SiO masers. This position is consistent 
with that estimated from the distribution of maser spots with an uncertainty better than 10 mas. 
We found there are several spoke-like maser features; most of these could be modeled by ballistic 
orbits, suggesting that most spoke-like maser features are on the far side of the star, are decelerating, 
pointing back to the central star. 

The kinematics of SiO maser spots show a slow expansion. After removing the effects of 
expansion, we derived an absolute proper motion of the central star. A large discrepancy of tens of 
mas between the positions of VY CMa determined by optical and radio measurements is found, and 
we suggest that optical position determinations are affected by the scattering from circumstellar 
dust and cannot be used to locate the star or its motion. 
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